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Abstract 

We present updated measurements of the branching fractions for the modes — > vr'^vr^, B^^ 
TT^vr", and B^ — > K^tt^. We also measure the time- integrated asymmetry C^o^o and the charge 
asymmetries .47r±7r0 and Ak±ttO- Based on a sample of approximately 227 million BB pairs collected 
by the BABAR detector at the PEP-II asymmetric-energy B Factory at SLAG, we measure 

fi(B° ^ 7r°7r°) = (1.17 ± 0.32 ± 0.10) x 10^^ , C^o^o = -0.12 ± 0.56 ± 0.06 , 

where the first errors are statistical and the second are systematic. The B^ — > vr'^vr'^ signal has a 
significance of 4.9 a including systematic uncertainties. We also measure 

B{B^ 7r=^7r°) = (5.8 ± 0.6 ± 0.4) x 10"^ , A^±^o = -0.01 ± 0.10 ± 0.02 , 
BiB^ K^ir°) = (12.0 ± 0.7 ± 0.6) x 10"^ , Ak±^o = 0.06 ± 0.06 ± 0.01 . 

Using related BABAR measurements and isospin relations we find an upper bound on the angle 
difference \6\ = \a — agfrl of 35° at the 90% C.L. 
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1 INTRODUCTION 



The study of B meson decays into charmless hadronic final states is important for the understanding 
of CP violation in the B system. In the Standard Model, CP violation arises from a single phase 
in the Cabibbo-Kobayashi-Maskawa quark-mixing matrix V^q/ pp. Measurements of the time- 
dependent CP-violating asymmetry in the B^ tt^tt^ decay mode by the BABAR and Belle 
collaborations [21 01 provide information on the angle a = arg [— l^d^b/^ud^ubl °^ Unitarity 
Triangle. However, in contrast to the theoretically clean determination of the angle (3 in B^ 
decays to charmonium final states IHEI, the extraction of a in B^ — > vr^vr" is complicated by the 
interference of tree and penguin amplitudes with different weak phases. The difference between 
the angles a and Oes, where a^s is derived from the measured time-dependent B'^ vr'^vr" CP 
asymmetry, may be evaluated using measurements of the isospin-related decays B^{B ) — > 
and B^ — > vr^vr*^ 0. The amplitudes A^^ for the B vrV-' decays satisfy the relation 

= -^A+- + (1) 
v2 

and a similar relation for the conjugate amplitudes A'^ , A^^ and ^4^'' = ^4"*^. 

For 5° TT^n^, the CP-related quantity that can be accessed experimentally is the time- 
integrated asymmetry C^o^o, defined as 

1 - IAqoI^ 
t + |Aoo| 



where lAi 



00 1 



Aqq /Aqq is the modulus of the ratio of B and B^ tt^tt^ decay amplitudes. 

C^o^o may deviate from zero if the tree and penguin amplitudes are of comparable size and have 
different weak and strong phases. For B^ modes, the CP-violating charge asymmetry is defined as 

_ \A\-- \A\- 

■^^^ - \A\^ + \Ar ^ ^ 

where A {A) is the B^ ) decay amplitude. In the Standard Model, the decay B^ — > vr^vr^ is 
governed by a pure tree amplitude; as a result no charge asymmetry is expected. Both the rate 
and asymmetry of the decay B^ —>■ K^tt^ may be used to extract useful constraints on penguin 
contributions to the B Kn amplitudes. 

In this paper, we report a measurement of the time-integrated CP asymmetry in the B^ tt^tt^ 
decay, and an updated measurement of the B branching fraction [3|. We also update the 

B^ TT^vr" and B^ — > K^ir^ branching fractions and charge asymmetries ISl. This study is based 
on (226.6 lb 2.5) x 10^ T(45) — > BB decays (on-resonance data), collected with the BABAR detector. 
We also use approximately 16 fb~^ of data recorded 40 MeV below the BB pair production threshold 
(off-resonance data). 



2 THE BABAR DETECTOR AND PARTICLE SELECTION 

BABAR is a solenoidal detector optimized for the asymmetric-energy beams at PEP-II and is 
described in detail in Ref. Charged particle (track) momenta are measured with a 5-layer 
double-sided silicon vertex tracker and a 40-layer drift chamber inside a 1.5-T superconducting 
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solenoidal magnet. Neutral cluster (photon) positions and energies are measured with an electro- 
magnetic calorimeter (EMC) consisting of 6580 CsI(Tl) crystals. The photon energy resolution 
is ge/E = {2.3/E(GeV)V4 i.gj 

%, and the angular resolution from the interaction point is 

ae = 3.9° /^E{ GeV). The photon energy scale is determined using symmetric 7r° — > 77 decays. 

Charged hadrons arc identified with a detector of internally reflected Cherenkov light (DIRC) and 
using ionization measurements in the tracking detectors. An instrumented magnetic-flux return de- 
tects neutral hadrons and identifies muons. High efficiency for recording BB events in which one B 
decays with low multiplicity is achieved with a two-level trigger with complementary tracking-based 
and calorimetry-based trigger decisions. 

Candidate tt^ mesons are reconstructed as pairs of photons, spatially separated in the EMC, 
with an invariant mass m^o satisfying 110 < 77^7^0 < 160 MeV/c^. The resolution is approximately 
8 MeV/c^ for high momentum tt'^'s. Photon candidates are required to be consistent with the 
expected lateral shower shape, not be matched to a track, and have a minimum energy of 30 MeV. 
To reduce the background from false Tr*^ candidates, the angle 9^ between the photon momentum 
vector in the vr*^ rest frame and the tt*^ momentum vector in the laboratory frame is required to 
satisfy |cos^-y| < 0.95. 

Candidate tracks are required to be within the tracking fiducial volume, originate from the 
interaction point, consist of at least 12 DCH hits, and be associated with at least 6 Cherenkov 
photons in the DIRC. 

3 ANALYSIS METHOD 

B meson candidates are reconstructed by combining a vr" with a pion or kaon (h^) or by combining 
two 7r° mesons. Two kinematic variables, used to isolate the B^ — > tt^tt^ and B^ —>■ h^n'^ signal 
events, take advantage of the kinematic constraints of B mesons produced at the T{AS). The first is 
the beam-energy-substituted mass uies = \j is/2 + Pj • t^bY /Ef — p^, where ^/s is the total e^e~ 
center-of-mass (CM) energy, (£'j,pj) is the four-momentum of the initial e^e~ system and is 
the B candidate momentum, both measured in the laboratory frame. The second variable is IS.E 
= Eb — \fs/2, where Eb is the B candidate energy in the CM frame. The resolution for signal 
is approximately 80 MeV for B^ 7r°7r°, and 40 MeV for B^ /t=^7r°. 

The primary source of background is e^e~ qq {q = u, d, s, c) events where a tt^ or from 
each quark randomly combines to mimic a B decay. The jet-like qq background is suppressed by 
requiring that the angle 6s between the sphericity axes of the B candidate and that of the remaining 
tracks and photons in the event, in the CM frame, satisfies |cos^s| < 0.7(0.8) for B'^ — ^ tt^tt^ 
{B^ h^TT^). The other sources of background are B to vector-pseudoscalar decays: B^ p^ir^ 
for the ^ tt^ttO mode; S° ^ p^-K^ , B^ p°7r=^, and B^ -> p=^7r° for the B"^ tt^tt^ mode; 
and 5° P^K^, B^ K*^n^, and 5° i^*°7r° for the B^ K^-k^ mode. In these decays the 
vector meson is polarized, so one of the pions is often produced almost at rest in the B rest frame, 
and the remaining decay products match the kinematics of a B^ 7r*^7r^ or B^ h^TT^ decay. 

For the B^ tt^tt^ analysis we restrict the w-es — ^E plane to the region with rriES > 5.2 
GeV/c^ and |Ai?| < 0.4 GeV. For the on-resonance sample we define the signal region as the 
band in the plane with |A£^| < 0.2 GeV and the sideband region as the rest of the plane with the 
exception of a region at negative AE also populated with B^ — > p^ir^ events. The entire plane 
for the off-resonance data and the sideband region for the on-resonance data are populated with qq 
background events, which are kept in the fit in order to constrain the corresponding background 
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parameters. /i^vr^ candidates are selected in the region with rngs > 5.22 GeV/c'^ and 

-0.11 <AE< 0.15 GeV. 

For — > vr'^vr'^ candidates, the other tracks and clusters in the event are used to determine 
whether the other B meson (Stag) decayed as a B° or B^ (flavor tag). We use a multivariate 
technique ^H] to determine the flavor of the Stag meson. Separate neural networks are trained to 
identify primary leptons, kaons, soft pions from D* decays, and high-momentum charged particles 
from B decays. Each event is assigned to one of several mutually exclusive tagging categories based 
on the estimated mistag probability and on the source of tagging information. 

The number of signal B candidates is determined in an extended unbinned maximum likelihood 
fit. The probability Vt {xj;di) for a signal or background hypothesis is the product of probability 
density functions (PDFs) for the variables Xj given the set of parameters di. The likelihood function 
is given by a product over the events and the M signal and background hypotheses: 



/ M \ N 

exp - 



n 

i=i / j=i 



■ M 
.i=l 



(4) 



For B^ vr'^vr'^ the coefficients Ni are given by Ni = ^(1 — SjAi)ni, where sj refers to the sign 
of the flavor tag of the other B in the event j and is zero for untagged events, rii and Ai are 
the number of events and raw asymmetry for B^ — > vr'^vr'^ signal, B^ — > p^ir^ background, and 
continuum background components. The averages of the B^ — > p^Tr^ branching fraction and 
asymmetry |lllll2j are used to fix the number of B~^ p'^n^ and B~ p~iT^ events to 32.1 it 5.6. 
The asymmetry A^o^o is related to the time-integrated CP asymmetry C^o^o by the relation 

Ao^o = (l-2x)(l-2a;)C^o^o, (5) 

where x = 0.186 it 0.004 13 is the time-integrated mixing probability, and uj is the mis-tagging 
probability. Asymmetries in mistag rates and efficiencies are taken into account. 

For B^ — > h^ir'^ the probability coefficients are iVj = ^(1 — qjAi)ni, where qj is the charge 
of the track h in the event j, and the fit parameters nj and Ai are the number of events and 
asymmetry for B^ — > tt^tt^ and B^ K^ir^ signal, continuum, and B background components. 
The numbers of B background events are fixed to the expected number of events using the current 
world averages for the pn |111 11211111 1151 [T^ and pK decays, which are 18 it 4 and 3ib 1 events, 
respectively. We estimate a contribution of 1 it 1 S — > K*Tr^ events to the background based on 
the upper limit on the branching fraction The uncertainties on these numbers are dominated 
by the uncertainty on selection efficiencies, because of the sensitivity to the tight requirement in 
AE. 

The variables xj used for B^ — > vr'^vr'' are rriEs, AE, and a Fisher discriminant J-. The Fisher 
discriminant is an optimized linear combination of J2i Pi ^^'^ J2i Pi cos^ 6i , where pi is the momentum 
and di is the angle with respect to the thrust axis of the B candidate, both in the CM frame, for all 
tracks and neutral clusters not used to reconstruct the B meson. For both the B^ — > vr'^Tr'' signal 
and the B^ — > p^n^ background the mES and AE variables are correlated and therefore a two- 
dimensional PDF from a smoothed, simulated distribution is used. For the continuum background, 
the TTT-ES distribution is modeled as a threshold function jJH], whose shape parameter ^ is free in the 
fit, and the AE distribution as a quadratic function whose parameters are free in the fit. The PDF 
for the J- variable is modeled as a parametric step function (PSF) 7. for all event components. A 
PSF is a binned distribution whose parameters are the heights of each bin. The PSF is normalized 
to one, so that the number of free parameters is equal to the number of bins minus one. The J- PSF 
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has ten bins chosen so that each bin contains 10% of the signal sample. For the — > vr'^vr^ signal 
and the p^vr" background the J- PSF parameters are correlated with the flavor tagging, and 

the PSF parameters are different for each tagging category. Simulation is used to determine the 
PSF distributions for both S° 7r°7r° and B^ p±7r°. A sample of B ^ D^*^mT (n = 1,2,3) 
events is used to verify that the simulation reproduces the distribution. For qq background, the 
J- PSF parameters are free parameters in the fit. 

The variables Xj used for B^ h^i:^ are tues, Ai?, the Cherenkov angle 9c of the track, 
and the Fisher discriminant J-. The PDF parameters for m-Es, AE, 6c, and J- for the background 
are determined using the data, while the PDFs for signal are found from a combination of simulated 
events and data. The mss and /S.E distributions for qq events are treated as in the B^ vr^vr'^ case, 
with parameters allowed to vary freely in the fit. For the signal, the tties and IS.E distributions 
are both modeled as a Gaussian distribution with a low-side power law tail whose parameters are 
determined from simulation. The mean of the signal rriES and /S.E distributions are determined 
from the fit to the B^ /i^vr" sample and their values used to tune the vr" energy scale in the 
B^ vr^vr^ analysis. The mean of /S.E for the B^ K^t:^ mode is a function of the kaon 
laboratory momentum, since a pion mass hypothesis is used. The distribution of J- is modeled as a 
bifurcated Gaussian for the signal, whose parameters are obtained from simulation, and as a double 
Gaussian for the continuum background, whose parameters are determined in the likelihood fit. 
The difference of the measured and expected values of 9c for the pion or kaon hypothesis, divided 
by the uncertainty on 9c-, is modeled as a double Gaussian function. A control sample of kaon and 
pion tracks, from the decay D*^ — > D^ir^ , — > i^^vr"*", is used to parameterize ao^ as a function 
of the charged track momentum. 

4 RESULTS AND SYSTEMATIC UNCERTAINTIES 

The result of the maximum likelihood fit for S° 7r°7r° is n{B^ 7r°7r°) = 61 ± 17 (see Table CJ). 
The statistical significance of the event yield is evaluated from the square root of the change in 
— 21n£ between the nominal fit and a separate fit in which the signal yield is fixed to zero, and is 
found to be 5.2cr (statistical errors only). The time-integrated asymmetry is C^a^o = —0.12 it 0.56. 
Distributions of likelihood fit variables for B^ — > vr'^vr'^ are shown in Fig. ^ The data shown are 
for events that have passed a probability ratio cut optimized to enhance the signal to background 
fraction. A validation is made by performing a simpler event-counting analysis, based on the 
number of B candidates satisfying tighter requirements. The result of the event-counting 

analysis is n{B^ — > vr'^vr'^) = 43ib26 and agrees well with the result from the maximum likelihood fit. 
This result is statistically consistent with our previously reported measurement J^. With changes 
in the analysis technique to measure the CP asymmetry, we now find 44 ± 13 events in the first 
123 million BB events, compared to 46 ± 13 events found in Ref. 0- The additional 104 million 
BB events dataset has a signal of 17 it 11 events. The signal rates in these two subsets agree at 
the 1.3(7 level. This result also includes an improved understanding of the vr*^ detection efficiency. 
Using a sample of vr^ mesons from decays we establish a vr efficiency correction of 

e° = 0.99 ± 0.03 to our GEANT simulation, compared to a correction of 0.88 ± 0.08 in Ref.0. 

For B^ h^-K^ the results are n{B^ 7r=^7r°) = 379±41 and n{B^ K^-k^) = 682±39. The 
charge asymmetries are At^±^q = —0.01 it 0.10 and ^^±^0 = 0.06 ± 0.06. Results are summarized 
in Tabled Distributions of likelihood fit variables for B^ h^ir^ are shown in Fig. [21 

Systematic uncertainties on the event yields and CP asymmetries are evaluated on data control 
samples, or by varying the fixed parameters and refitting the data. Tables |21 and El summarize 
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the dominant contributions to the systematic uncertainties for the — > vr'^vr'^ and /i± IT 

branching fractions and asymmetries. For B^ — > vr'^vr'', the dominant systematics arise from the 
uncertainty on the AE resolution and the efficiency of the vr^ reconstruction. We reevaluate the 
significance of the B vrOvrO si gnal from zero, including all systematic effects in the direction 
which reduces the signal, and find a significance of 4.9o". 

For both B^ tt^tt^ and B^ K^n^, the dominant systematic uncertainties arise from the 
PDF parameters for signal, selection efficiencies, and the AE resolution. Additional systematic 
uncertainties arise from the EMC energy scale, selection efficiencies, and particle identification. The 
systematic uncertainty on the charge asymmetries is dominated by the upper limit of the charge 
bias in the detector (1.0%) ^j, since most selection and PDF systematics cancel in the asymmetry. 

We use the isospin relations, Eq.^ to extract information on the angle difference 6 = a — Oeff [S]; 
in conjunction with the BABAR measurements of 0^^+^- = (—0.09 it 0.15 it 0.04) j2j, the branching 
fraction i3(5° vr+vr-) = (4.7 ±0.6 ±0.2) x 10"^ [13], the S° 7r°7r° and B^ tt^tt^ decay rates 
and the 6*^0,^0 asymmetry described in this paper. We scan over all values of |<5| and calculate a 
for the five amplitudes {A'^ , A^ , A'^^, A and A'^^), given these five measurements and the two 
constraints (Eq. n^-nd its conjugate) for each value of \5\. The is converted into a confidence 
level, as shown in Fig. |31 from which we derive an upper bound on \5\ of 35° at the 90% C.L. 

5 CONCLUSION 

We observe 61 ± 17 ± 5 B^ tt^tt^ events with a significance of 4.9 a including systematic uncer- 
tainties. This corresponds to a branching fraction B{B^ — > vr'^vr'^) = (1.17 ± 0.32 ± 0.10) x 10~^, 
where the first error is statistical and the second is systematic. We measure the time-integrated 
asymmetry Ct^o^^q = —0.12 ± 0.56 ± 0.06. We report branching fractions B{B^ — > tt^tt^) = 
(5.8 ± 0.6 ± 0.4) X 10"*^ and B{B^ K^tt^) = (12.0 ± 0.7 ± 0.6) x lO"*^. The charge asym- 
metries are A^±t^o = —0.01 ± 0.10 ± 0.02 and ^^^±7^0 = 0.06 ± 0.06 ± 0.01; we find no evidence for 
CP violation. These results are consistent with our previous results for these decays [71|5]. We find 
an upper bound on the angle difference 5 = a — CKeg oi\5\ < 35° at the 90% C.L. 
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Table 1: The results for the modes B'^ —>■ tt^tt^ and h^n'^ are summarized. For each mode, 

the sample size N, number of signal events Ng, total detection cfRcicncy e, branching fraction B, 
asymmetry A or C^Oj^o, and 90% confidence interval for the asymmetry are shown. For C^o^o the 
confidence interval is normalized to the physical region [—1, 1]. Errors are statistical and systematic 
respectively, with the exception of e whose error is purely systematic. 



Mode 


N 


Ns 


£(%) 


i3(10-6) 


Asymmetry 


(90%C.L.) 


B^^ 




8153 


61 ±17 


23.5 ± 1.4 


1.17 ±0.32 ±0.10 


-0.12 ±0.56 ±0.06 


[-0.88, 0.64] 


B^^ 


7r±7rO 


29950 


379 ± 41 


28.7 ± 1.1 


5.8 ±0.6 ±0.4 


-0.01 ±0.10 ±0.02 


[-0.19,0.21] 


B^^ 




13165 


682 ± 39 


25.0 ± 1.0 


12.0 ±0.7 ±0.6 


0.06 ± 0.06 ± 0.01 


[-0.06,0.18] 



Table 2: Systematic uncertainties in the determination of the B^ — > branching fraction (left), 
and the C^o^o asymmetry (right). In order of decreasing importance the branching fraction is 
affected by the uncertainties in the efficiency (3% for each tt"), in the difference between simulation 
and data of the AE resolution, in the branching fraction B"^ P^7i"°, and in the number of 
produced BB pairs. The systematic uncertainty on C^o^o is dominated by the uncertainty on the 
B background asymmetry. 



Source; 






S(Mirce 


^(CV-'tt'') 


vr'' efficiency 


±6% 


-6% 


B background asymmetry 


±0.05 


AE resolution 


±7.2% 


-1.0% 


tagging efficiency 


±0.03 


B{B pn) 


+3.4% 


-3.1% 


simulation 


±0.02 


mean of AE and uies 


±1.5% 


-1.6% 


B{B^ Ptt) 


±0.01 


luminosity 


±1.1% 


-1.1% 







Table 3: Dominant systematic uncertainties for B^ — > tt^tt'^ and B^ —>■ K^i:^, listed in order of 
decreasing importance as percent changes in the branching fractions B (left), and absolute changes 
in the asymmetries A^±j^o, AK±n° (right)- 



Source 


AB(7r±7r°) 


AH(K±7r°) 


Source 


A(A±^o) 


A{AK±no) 


AE simulation 


±4.2% 


±3.3% 


detector bias 


±0.010 


±0.010 


7r° efficiency 


±3.1% 


±3.1% 


AE simulation 


±0.011 


±0.007 


B background 


±2.2% 


±0.2% 


B background 


±0.008 


±0.001 


luminosity 


±1.1% 


±1.1% 


identification 


±0.004 


±0.003 


identification 


±0.5% 


±0.6% 
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Figure 1: The distributions of (a) mES? (b) AE, and (c) Fisher discriminant for — > vr^vr'^ 
candidates in the signal data sample that satisfy an optimized requirement on the signal probability, 
based on all variables except the one being plotted. The projections contain 25%, 45% and 68% of 
the signal, 14%, 31% and 17% of the pir^ background, and 2.2%, 1.3% and 4.4% of the continuum 
background, for m-Es, AE, and J- respectively. The PDF projections are shown as a dashed line for 
qq background, a dotted line for — > /O^vr*^, and a dashed-dotted line for B^ — > vr'^Tr'^ signal. The 
solid line shows the sum of all PDF projections. The PDF projections are scaled by the expected 
fraction of events passing the probability-ratio requirement. 
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Figure 2: The distributions and PDF projections of mss (top), IS.E (middle) and J- discriminant 
(bottom), for tt^-k^ and — > K^tt^ candidates in the data sample. Figures on the left are 

for background-subtracted signal events, and background distributions are on the right. For rriEs 
and J- discriminant, the B^ — > vr^vr'^ and B^ — > K^tt^ distributions are combined, while for signal 
the B^ — > tt^tt'' (open circles and dotted dashed curve) and B^ — > K^tt'^ (solid circles and 
solid curve) distributions are shown separately. The background is subtracted using the method 
described in reference [23|. The method uses all variables except the one being plotted. 
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Figure 3: Constraints on the angle 5 = a — aeff, expressed as a Confidence Level as a function of 
\6\. The constraint is evaluated using the isospin relations (Eq. ^) and the BABAR measurements 
of the B vrvr branching fractions and the CP violation parameters Cj^o^o and C^+^- . We find an 
upper bound on |(5| of 35° at the 90% C.L. 
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